Isospin mixing in the nucleon and 4 He and the nucleon strange electric form factor 
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In order to isolate the contribution of the nucleon strange electric form factor to the parity- 
violating asymmetry measured in 4 He(e, e') 4 He experiments, it is crucial to have a reliable estimate 
of the magnitude of isospin-symmetry-breaking (ISB) corrections in both the nucleon and 4 He. We 
examine this issue in the present letter. Isospin admixtures in the nucleon are determined in chiral 
perturbation theory, while those in 4 He are derived from nuclear interactions, including explicit 
ISB terms. A careful analysis of the model dependence in the resulting predictions for the nucleon 
and nuclear ISB contributions to the asymmetry is carried out. We conclude that, at the low 
momentum transfers of interest in recent measurements reported by the HAPPEX collaboration at 
Jefferson Lab, these contributions are of comparable magnitude to those associated with strangeness 
components in the nucleon electric form factor. 
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One of the challenges of modern hadronic physics is to 
determine, at a quantitative level, the role that quark- 
antiquark pairs, and in particular ss pairs, play in the 
structure of the nucleon. Parity-violating (PV) electron 
scattering from nucleons and nuclei offers the opportu- 
nity to investigate this issue experimentally. The PV 
asymmetry (Apy) arises from interference between the 
amplitudes due to exchange of photons and Z-bosons, 
which couple respectively to the electromagnetic (EM) 
and weak neutral (NC) currents. These currents involve 
different combinations of quark flavors, and therefore 
measurements of Apy, in combination with electromag- 
netic form factor data for the nucleon, allow one to iso- 
late, in principle, the electric and magnetic form factors 
G 8 E and G S M , associated with the strange-quark content 
of the nucleon. 

Experimental determinations of these form factors 
have been reported recently by the Jefferson Lab 
HAPPEX [l| and GO Collaborations, Mainz A4 
Collaboration [|, and MIT-Bates SAMPLE Collabora- 
tion These experiments have scattered polarized 
electrons from either unpolarized protons at forward an- 
gles 0, 0, IH or unpolarized protons and deuterons at 
backward angles Q. The resulting PV asymmetries are 
sensitive to different linear combinations of G S E and G S M 
as well as the nucleon axial- vector form factor G Z A . How- 
ever, no robust evidence has emerged so far for the pres- 
ence of strange-quark effects in the nucleon. 

Last year, the HAPPEX Collaboration 0, @] at Jeffer- 
son Lab reported on measurements of the PV asymmetry 
in elastic electron scattering from 4 He at four-momentum 
transfers of 0.091 (GeV/c) 2 and 0.077 (GeV/c) 2 . Because 
of the J n =0 + spin-parity assignments of this nucleus, 
transitions induced by magnetic and axial- vector currents 



are forbidden, and therefore these measurements can lead 
to a direct determination of the strangeness electric form 
factor G S E 7, 8], provided that isospin symmetry breaking 
(ISB) effects in both the nucleon and 4 He, and relativis- 
tic and meson-exchange (collectively denoted with MEC) 
contributions to the nuclear EM and weak vector charge 
operators, are negligible. A realistic calculation of these 
latter contributions [8j found that they are in fact tiny at 
low momentum transfers. The goal of the present letter 
is to provide a quantitative estimate of ISB corrections 
to the PV asymmetry. 

In the following analysis, we only need to consider the 
time components of the EM current and vector part of 
the weak NC current — the weak vector charge referred 
to above Q . We account for isospin symmetry breaking 
in both the nucleon and a-particle. We first discuss it in 
the nucleon. 

Ignoring radiative corrections, the EM and weak vector 
charge operators can be decomposed as 
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where j ^ ' and j'W are respectively the isoscalar and 
isovector components of the EM charge operators, j'W 
is the (isoscalar) component due to strange-quark contri- 
butions, and s^y — sin 2 6w contains the Weinberg mixing 
angle. In a notation similar to that adopted by the au- 
thors of Ref. [9( , we introduce form factors corresponding 
to the following matrix elements of and between 
proton (p) and neutron (n) states: 
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(p\j {1) \p) - G^(g 2 ) + 4(g 2 ) , (5) 

(n\j^\n) - -G^(g 2 ) + 4(Q 2 ) , (6) 

where the arrow indicates that only leading contributions 
are listed in the non-relativistic limit of these matrix el- 
ements. While higher order corrections associated with 
the Darwin- Foldy and spin-orbit terms are not displayed 
explicitly in the equations above, they are in fact re- 
tained in the calculations discussed later in the present 
work. The form factors G E (Q 2 ) and G E (Q 2 ) parameter- 
ize ISB effects in the nucleon states. We also introduce 
the strange form factor via 



(p\j {s) \p)=(n\j {s) \n)^G%(Q 2 ) 



(7) 



where here ISB terms in the p,n states are neglected. 
Contributions from sea quarks heavier than strange are 
also ignored. 

In terms of the experimental proton and neutron 
electric form factors, derived from the matrix elements 
Mm» - G E (Q 2 ) and (n\j£=°\n) - G E (Q 2 ), we ob- 



tain: 
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G E = (G P E - G E )/2 - G E 



(8) 
(9) 



where the Q 2 dependence in these and the following two 
equations is understood. In the limit in which the p,n 
states form an isospin doublet, the form factors G E and 
G E vanish, and G E and G\ reduce to the standard 
isoscalar and isovector combinations of the proton and 
neutron electric form factors. The proton and neutron 
vector NC form factors follow from Eq. (J2]), i.e. 

G P E Z = (1 - 4s 2 w )G E -G% + 2(G E ~ G%) - G% , (10) 

G" E Z = (1 - ls 2 w )Gl -G E + 2(4 + G%) - G% . (11) 

We now turn to the nuclear charge operator. At low 
momentum transfer, it is simply given by 
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where Z is the number of protons, A — Z the number of 
neutrons, and for elastic scattering from a nuclear target 
of mass niA the squared four-momentum transfer is taken 
as Q 2 — 2 niA ( \/ q 2 + m 2 A — ttla ) , with q being the three- 
momentum transfer, and q = |q|. An equation similar to 
Eq. (fT2| holds for the weak vector charge operator, but 
with G P E and G E being replaced respectively by G P E and 
G n E . It is also convenient to define the charge operators: 
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p^(q) =pW(q)+pW(q), 
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where again the Q 2 dependence of the nucleon form fac- 
tors has been suppressed here and in the following for 
brevity. The relations above lead to the definition of the 
following nuclear form factors: 

( A Re\p^(d)fRe)/Z = F^( q ) , a = EM, 0,1, (17) 

having the normalizations i^( EM ) (0)=F(°> (0)=1 and 
FW(0)=0. The form factor F^(q) is very small be- 
cause He is predominantly an isoscalar state. Thus, ig- 
noring second order terms like G FW(«), we obtain for 
the PV asymmetry measured in (e, e') elastic scattering 
from 4 He: 
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where G M is the Fermi constant as determined from muon 
decays, and here s\y is taken to incorporate radiative 
corrections. The terms G E and F^(q)/F^(q) are the 
contributions to Apy, associated with the violation of 
isospin symmetry at the nucleon and nuclear level, re- 
spectively. 

The most accurate measurement of the PV asymmetry, 
recently reported in Ref. 6] at Q 2 =0.077 (GeV/c) 2 , gives 
A PV = [+6.40 ± 0.23 (stat) ± 0.12 (syst)] ppm, from 
which, after inserting the values for 0^=1.16637 x 10 -5 
GeV" 2 , a=l/137.036, and s^=0.2286 (including its ra- 
diative corrections |7J) in Eq. (JT5J) , one obtains 



2G' 



G s , 
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0.010 ± 0.038 (19) 



at Q 2 = 0.077 (GcV/c) 2 . This result is consistent with 
zero. In the following, we discuss the estimates for the 
ISB corrections first in the nucleon and then in 4 He, re- 
spectively G E {Q 2 ) and at Q 2 =0.077 (GeV/c) 2 
(corresponding to q=lA fin -1 ). 

For G E (Q 2 ) we use the estimate obtained in Ref. [§] 
adapted to our conventions, combining a leading-order 
calculation in chiral perturbation theory with estimates 
for low-energy constants using resonance saturation. Col- 
lecting the various pieces, we find 
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FIG. 1: The isospin-violating nucleon form factor G^ E (Q 2 ). 
The band comprises a range of values for various vector-meson 
coupling constants, as well as an estimate of higher-order chi- 
ral corrections. For details, see Ref. 
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where the loop functions £, 7 / 3 are given explicitly in 
Ref. along with the precise definitions of the vari- 
ous coupling constants. The chiral loop contributions 
in Eq. (|20|) scale with the neutron-proton mass differ- 
ence Am, while the resonance part is proportional to the 
p-u mixing angle 9 paJ . We refer to Ref. [§] for a de- 
tailed discussion of the range of numerical values for the 
vector meson coupling constants and only show the re- 
sulting band for G E {Q 2 ) in Fig. [TJ At the specific kinc- 
matical point of interest Q 2 =0.077 (GeV/c) 2 , we find 



G E (Q 



0.0017 ± 0.0006, and with G^(Q 2 )=0.799 



and G|(Q 2 )=0.027 0, we obtain 



2G\ 



(G P E + G%)/2 



= 0.008 ±0.003 



(21) 



at Q 2 = 0.077 (GeV/c) 2 . 

We now turn to the nuclear ISB corrections. An ap- 
proximate calculation of the ratio F^(g)/F^°\q) was 
carried out more than a decade ago [ll|, by i) tak- 
ing into account only the isospin admixtures induced by 
the Coulomb interaction, ii) constructing a T=l J 7r =0 + 
breathing mode excitation based on a plausible ansatz, 
and iii) generating the relevant T=\ component in the 
4 Hc ground state in first order perturbation theory. The 
calculated value was found to be rather small, and it pro- 
duced a less than 1% correction with respect to the 4.s 2 v 
term in Eq. (02) at low Q 2 . 

Since that pioneering study, significant progress has 
occurred on several fronts. First, there now exist a num- 
ber of accurate models of nucleon-nucleon (NN) poten- 



tials [H, El El El Ell which include explicit ISB in- 
duced by both the strong and electromagnetic interac- 
tions. These ISB terms have been constrained by fitting 
pp and np elastic scattering data. It is now an established 
fact that a realistic study of 4 He, and in fact light nu- 
clei [TtJ , requires the inclusion of three- nucleon (NNN) 
potentials in the Hamiltonian. While these are still not 
well known, the models most commonly used in the lit- 
EI El 0| do not contain ISB terms. The 
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strength of the latter, however, is expected to be tiny. 

Second, several accurate methods have been developed 
to compute 4 He wave functions starting from a given re- 



alistic nuclear Hamiltonian |2l|. In these calculations, 
T > components are generated non-perturbatively. 
The T=l percentage in the 4 He wave function is typi- 
cally found to be of the order of 0.001 %. 

In this paper, we use the Hyperspherical Harmonic 
(HH) expansion method to compute the 4 He wave func- 



tion [22], [23, |24|. In order to have an estimate of the 



model dependence, we consider a variety of Hamilto- 
nian models, including: i) the Argonne t> ls NN poten- 
tial 0JAV18); ii) the AV18 plus Urbana-IX NNN po- 
tential (AV18/UIX); iii) the CD Bonn Q NN plus 
Urbana-IXb NNN potentials (CDB/UIXb); and iv) the 
chiral N3LO NN potential (N3LO). The Urbana 
UlXb NNN potential is a slightly modified version of 
the Urbana UIX (in the UlXb, the parameter Uq of the 
central repulsive term has been reduced by the factor 
0.812), designed to reproduce, when used in combination 
with the CD Bonn potential, the experimental binding 
energy of 3 H. The binding energies B and Pt=i per- 
cent probabilities obtained with the AV18, AV18/UIX, 
CDB/UIXb, and N3LO are respectively £=(24.21, 28.47, 
28.30, 25.38) MeV (to be compared with an experi- 
mental value of 28.30 MeV) and P T =i = (0.0028, 0.0025, 
0.0020, 0.0035). These results are in agreement with 
those obtained with other methods (for a comparison, 
see Ref. j^). 

The form factors F^°'(q) and F^'(q), defined in 
Eq. (JTTJ) and calculated with the AV18/UIX Hamiltonian 
model, are displayed in Fig. [2] The dashed (solid) curves 
represent the results of calculations including the one- 
body (one-body plus MEC) EM charge operators (note 
that ISB corrections in the nucleon form factors enter- 
ing the two-body EM charge operators, listed explicitly 
in Ref. Q, are neglected). Similar results (not shown 
in Fig. [2] to reduce clutter) are obtained with the other 
Hamiltonian models. In particular, the model depen- 
dence in the calculated F^ (q) form factor is found to be 
weak, although the change of sign in the predictions cor- 
responding to the N3LO model occurs at a slightly lower 
value of momentum transfer than in those corresponding 
to the other models, which are in excellent agreement 
with the experimental data from Refs. [2^]. From the 
figure it is evident that for q < 1.5 fm _1 , the effect of 
MEC in both F^ (q) and F^ (q) is negligible. 
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FIG. 2: The F (0) (q) and F (1) (q) form factors for the 
AV18/UIX Hamiltonian model. The F m (q) is compared 
with the experimental 4 He charge form factor [25l |. The ratio 
\F (1> (q) / F (0) (q)\ (all calculations include MEC) is shown in 
the inset for the four Hamiltonian models considered in this 
paper. 



In the inset of Fig. [2j we show the model depen- 
dence of the ratio \F^\q) / F^°\q)\ (all calculations in- 
clude MEC). The various Hamiltonian models give pre- 
dictions quite close to each other, although the value for 
the N3LO is somewhat larger than for the other models, 
reflecting the larger percentage of T=l admixtures in the 
4 He ground state, predicted by the N3LO potential. The 
calculated ratios F« (q)/F^ (q) at Q 2 =0.077 (GeV/c) 2 
are of the order of —0.002. The inclusion of NNN poten- 
tials tends to reduce the magnitude of F^/F^°\ while 
ignoring MEC contributions, at this value of Q 2 , would 
lead, at the most, to 1.5% decrease of this magnitude. 

Note that the value estimated in Ref. 
\F^/F^\ « 0.0014 at Q 2 =0.077 (GeV/c) 2 , 



11] 



although 

it was computed in first order perturbation theory by 
only keeping the ISB corrections due to the Coulomb po- 
tential. However, the latter only account for roughly 50 
% of the Pt=i probability in the 4 He ground state [23j |. 
and, assuming the ratio above to scale with \JPt=\, one 
would have expected a smaller value for it than actually 
obtained (» 0.0014) in Ref. [Il|. 

Therefore, at Q 2 =0.077 (GeV/c) 2 , both contributions 
F (D/ F (o) and G t E 

are found of the same order of magni- 



tude as the central value of T in Eq. (I19|) . Using in this 
equation the value F^/F^ » -0.00157 obtained with 
the Hamiltonian models including NNN potentials, and 
the chiral result for G E = -0.0017 ± 0.0006, one would 
obtain G% [Q 2 = 0.077 (GeV/c) 2 ] = -0.001 ± 0.016 thus 
suggesting that the value of T is almost entirely due to 
isospin admixtures. Of course, the experimental error on 
r is still too large to allow us to draw a more definite con- 
clusion. A recent estimate of G S E using lattice QCD input 
obtains [26| G%[0.1 (GeV/c) 2 ] = +0.001 ± 0.004± 0.003. 



An increase of one order of magnitude in the experimen- 
tal accuracy would be necessary in order to be sensitive 
to G S E at low values of Q 2 . Indeed, if the lattice QCD 
prediction above is confirmed, the present data would 
suggest that the leading correction to the PV asymme- 
try is from isospin admixtures in the nucleon and/or 4 He. 
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